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ABSTRACT 

Dark energy may be the first sign of new fundamental pliysics in the Universe, taking either a 
physical form or revealing a correction to Einsteinian gravity. Weak gravitational lensing and 
galaxy peculiar velocities provide complementary probes of General Relativity, and in com- 
bination allow us to test modified theories of gravity in a unique way. We perform such an 
analysis by combining measurements of cosmic shear tomography from the Canada-France 
Hawaii Telescope Lensing Survey (CFHTLenS) with the growth of structure from the Wig- 
gleZ Dark Energy Survey and the Six-degree-Field Galaxy Survey (6dFGS), producing the 
strongest existing joint constraints on the metric potentials that describe general theories of 
gravity. For scale-independent modifications to the metric potentials which evolve linearly 
with the effective dark energy density, we find present-day cosmological deviations in the 
Newtonian potential and curvature potential from the prediction of General Relativity to be 
A*/* = 0.05 ± 0.25 and A$/$ = -0.05 ± 0.3 respectively (68 per cent CL). 

Key words: cosmology: observations - gravitational lensing 
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1 INTRODUCTION 

Einstein's original formulation of General Relativity (GR) provides 
a remarkably precise prescription for the motions of particles in 
our solar system jEinstein|1916[ l. It has survived close to a century 
of experimental scrutiny, beginning with the precession of Mer- 
cury's orbit and the gravitational deflection of starlight, progressing 
to ground-based tests with atomic clocks. Recent spaceborne tests 
now confirm certain predictions of GR to better than one part in 
ten thousand ( [Bertotti et al.|2003t . However when we consider the 
behaviour of gravity on cosmological scales, which are over four- 
teen orders of magnitude greater than the inter-planetary distances, 
there is much greater scope for uncertainty. Observational evidence 
of an accelerating universe has called Einstein's laws into question. 
The underlying cause for this "dark energy" phenomenon may lie 
within our current understanding of gravitational physics, the lead- 
ing candidate being the cosmological constant. Another possibility 
is that a new regime of gravitational physics has been exposed. It 
is this prospect which has led to a recent flurry of activity in con- 



straining the cosmological nature of gravity (|Daniel et al. 
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Much like in the solar system, we are able to study the cos- 
mological trajectories of both relativistic and non-relativistic parti- 
cles, at least in a statistical sense. The peculiar motions of galaxies 
falling towards overdense regions generate an illusory anisotropy 
in their clustering pattern ( |Kaiser|1987| l, known as redshift space 
distortions (RSD). This effect has been observed with heightened 
precision over the past decade ("Peacock et al. 2001 ; Hawkins et al. 
2003||Guzzo et al.|2008||Blake et al.,201 l,.Beutler et al..2012„Reid, 
et gJ.|2012) permitting inferences of the rate at which cosmological 
structure forms. The peculiar velocities of galaxies can also be ob- 
served directly, and the inferred measurements of the growth rate 
are consistent with those from RSDs (see |Hudson & Tumbull|2012| 
and references therein). 

For gravitational experiments, relativistic test particles are of 
particular interest since their motions are not only sensitive to the 
effects of time dilation, but also to the curvature of space. In the 
classical test of GR, the deflection of starlight was observed during 
a solar eclipse jPyson et al.|192 0V In our experiment, distant galax- 
ies replace nearby stars as the source of photons, our relativistic 
test particle. When imaging these distant galaxies, instead of mea- 
suring their gravitational deflection angles, which are unknown, we 
can study their shapes which appear correlated as a result of the 
gravitational deflection of light i jBartelmann & Schneider|2001[ l. 

One important caveat is that we are required to convert an- 
gles and redshifts into physical distances when interpreting obser- 
vations of redshift space distortions and gravitational lensing. We 
therefore require additional information on the geometry of the uni- 
verse. Geometric measurements of dark energy such as supernovae 
and baryon acoustic oscillations, conventionally used to determine 
the equation of state for dark energy, are crucial for this pur- 

pose. However their measurements alone cannot distinguish mod- 
ified gravity models, which are in general capable of reproducing 
any given wCDM expansion history. The Cosmic Microwave Back- 
ground (CMB) is also used to further break degeneracies, and pro- 
vides a high redshift anchor for the amplitude of density perturba- 
tions (see e.g. Kom atsu et al.|2011| l. 

Gravitational lensing is not unique in its ability to probe both 
metric potentials. CMB photons gain or lose energy as the poten- 



tials evolve, a phenomenon known as the Integrated Sachs Wolfe 
effect (ISW), leaving a significant anisotropy on the largest angu- 
lar scales. In the context of a flat ACDM background, the expected 
signal strength from GR is modest, while many alternative theories 
of gravity can readily generate very large ISW signals, something 
which is not seen in the data ( |Lombriser et al.|2009| |Zuntz et al.| 
[2012) . 

While we currently lack a specific model which can compete 
with GR in terms of simplicity and physical motivation, there ex- 
ist broad families of models such as f{R) which provide useful 
test cases. For a recent review of the full menagerie of models, 
see |Clifton et al.| ( |2012| >. In this work we do not aim to constrain 
the parameter space of a particular family of models, or explic- 
itly perform Bayesian model selection, but instead simply address 
the question of whether our data appear consistent with the predic- 
tions of GR. In addition to the cosmic shear data from the Canada- 
France-Hawaii Telescope Lensing Survey (CFHTLenS), we make 
use of redshift space distortions from WiggleZ and the Six-degree- 
Field Galaxy Survey (6dFGS). 

In *j2]we outline our choice of parameterisation for departures 
from GR. Our datasets are summarised in 33] and Sj4] followed by 
our methods and main results in !|5] and Sj6] A comparison with 
alternative parameterisation schemes is presented in !|7] and we 
briefly review the status of theoretical models of modified gravity 
in ^ with concluding remarks in ^ 



2 PARAMETERISING THE MODIFICATION 

Given the lack of a compelling model to rival GR, we choose to 
parameterise deviations from GR in a phenomenological manner. 
In effect, we simply wish to ask: is the strength of gravity the same 
on cosmological scales as it is here on Earth? If not, this may mod- 
ify the motions of both relativistic and non-relativistic particles, but 
not necessarily in the same manner. Figure[T|illustrates why we are 
seeking a modified gravity signal on cosmological scales; it is the 
regime in which the classical Newtonian attraction is overwhelmed 
by the repulsive force associated with dark energy. 

The perturbed Friedmann-Robertson- Walker metric may be 
expressed in terms of the scale factor a(t), the Newtonian poten- 
tial ^, and the curvature potential $ 

ds^ = (1 + 2*) dt^ - a{t) (1 - 2$) dx^ . (1) 

In GR each of the two gravitational potentials and <J>, which 
carry implicit spatial and temporal dependencies, may be deter- 
mined from the distribution of matter in the Universe. In Fourier 
space this is given by 

k'^^GR = -^T^Ga?p5, (2) 

where we have introduced the wavenumber k, the mean cosmic 
density p, and the fractional density perturbation 5 = p/p ~ 1. 
In the absence of anisotropic stress the two potentials are identi- 
cal, '^{x,t) — ^{x,t). Relativistic particles collect equal contri- 
butions from the two potentials, since they traverse equal quantities 
of space and time, such that 

fc' (*GR + $gr) = -STvGa^pS . (3) 

We wish to investigate deviations from GR with a general parame- 
terisation. Hereafter, each of the two probes, the Newtonian poten- 
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tial ^ experienced by non-relativistic particles and the lensing po- 
tential ($ + vp) experienced by relativistic particles, are now mod- 
ulated by the parameters S and ^ 

^(fc.a) = [l + M(fc,a)]*GR(fc,ffl), (4) 



[*(fe, a) + ^(fc, a)] = [1 + E(fc, a)\ [^GYi{k, a) + $GR(fc, a)] , 

(5) 

where we have adopted notation similar to that of 'Amendola et al. 
1 2008) , as used more recently by [Zhao et al . (2010) and Song] 
|et al.| ( (201 1| |. This parameterisation has the advantage of separat- 
ing the modified behaviour of non-relativistic particles, as dictated 
by ii{k, a), from modifications to the deflection of light, as given 
by E(fe, a). By being able to reproduce a wide range of observa- 
tional outcomes, in terms of the two-point weak lensing correla- 
tion functions and the growth of large scale structure, we ensure a 
broad sensitivity to different types of deviations from GR. There 
is some flexibility in how we choose to parameterise the scale and 
time dependence of these two parameters, much like the dark en- 
ergy equation of state w[z). Previous works have often chosen a 
scale independent model with a parameterised time variation of 
these functions such that E(a) = Eso" and /i(a) = ^sO," . How- 
ever given that our primary motivation for modified gravity is to 
seek an explanation of the dark energy phenomenon, we model the 
time-evolution to scale in proportion with the effective dark energy 
density implied by the background dynamics, such that 

E(a) = Eo , /i(a) = , (6) 

where we have defined SIa = O.A{a — 1); this normalisation was 
chosen such that fio and Eo reflect the present day values of ^{a) 
and E(a) respectively. Note that for the case of GR, Eq = /^o = 0. 
One advantage of this parameterisation is that it permits a trivial 
mapping between /io and the popular growth index 7 jWang &| 
|Steinhardt|1998|[ljn der 2005). It also reduces the total number of 
degrees of freedom, thereby allowing us to maintain a manageable 
error budget. The disadvantage is that we become less sensitive to 
deviations from GR whose time evolution differs significantly from 
the adopted model. The lack of scale dependence in this model may 
seem questionable given the stringent Solar System tests which 
must be satisfied. However the length scales probed by our cosmo- 
logical datasets are over fourteen orders of magnitude greater, leav- 
ing ample opportunity for a transition at scales smaller than those 
studied here. Furthermore, any significant departure from scale- 
independent growth is more readily apparent from studies of the 
shape of the galaxy power spectrum. 

With our modified potentials in equations (|4]|5| the growth of 
linear density perturbations is now given by 



50. 




Distance (lO^in) 



Figure 1. The .shaded triangle represents the vast regime over which the 
gravitational acceleration of a test particle in the proximity of a mass M is 
thought to be well described by the Newtonian prescription a = GMr~^ . 
While torsion balance (Laboratory) and space-based (Solar System) tests of 
gravity offer superior accuracy, cosmological methods push the boundaries 
of our understanding of gravity. For clarity we truncate the mass axis at the 
Planck mass (red edge), while the upper bound (black edge) denotes the 
Schwarzschild radius. The grey slope illustrates the amplitude of the nega- 
tive contribution arising from dark energy, which scales in proportion to r, 
and surpasses the conventional Newtonian force at cosmological scales. 

where the terms not present in the equations' conventional form are 
highlighted in bold. Here x denotes the radial coordinate distance, 
while Jk is the comoving angular diameter distance. Note also that 
Pk {£) possesses an implicit dependence on fi via the time-evolution 
of the Ps term. The lensing efficiency gi is determined by the ra- 
dial distribution of source galaxies ni{x) and the comoving angular 
diameter distance /k (x) 

5.(X) = r dx'n.(x')^^44^ • (9) 
Jx Jk(x) 

Note that our parameters fio and Eq are only modulating the 
gravitational potentials in equation (jTJ, but not the evolution of the 
expansion a{t). The standard cosmological model, in the form of a 
flat ACDM universe, is known to provide a good description of the 
cosmic geometry, therefore for most of this paper we work on the 
assumption of this global expansion, and focus on exploring the 
less well determined behaviour of the gravitational perturbations. 
We shall also explore the consequences of adding further degrees 
of freedom to the expansion history, in the form of an effective dark 
energy equation of state w and the global curvature Qk- 



while the cosmic shear power spectrum P^'^ {£) correlating redshift 
bins i and j takes the form 

X [1 + E(x)]^dx, 

(8) 



3 GEOMETRIC DATA 

When considering the prospect of modified gravity, we must re- 
lax the assumption that light deflection and structure formation fol- 
low the predictions of GR. This significantly exacerbates the task 
of constraining all of our unknown cosmological parameters. For 
example, without knowing how mass bends light on cosmological 
scales, we cannot hope to use gravitational lensing to measure the 
geometry of the Universe. We therefore require auxiliary data to as- 
sist in measuring the conventional cosmological parameters, those 
which govern the cosmic expansion history and shape of the matter 
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power spectrum, before we can attempt to tackle the gravitational 
parameters jio and Eq. 



3.1 Hubble Parameter 

The current rate of cosmological expansion has long been a source 
of great uncertainty, largely due to the inherent difficulty in per- 
forming distance measurements. Considerable improvements have 
been made in the past decade, and we make use of a recent re- 
sult obtained from combining observations of local supemovae, 
Cepheid variables, and the megamaser at the centre of NGC 4258 
( |Riess et al.||201l"| >. Where stated, we include a Gaussian prior 



Riess et al. 



( |201l| >. In de 



Hq = 73.8 ± 2.4kms-^ Mpc"' from 
Sitter space this would correspond to the distsince between any two 
comoving points doubling every thirteen billion years. 

Since it is common practice to assume OR in most astronomi- 
cal and cosmological studies, extra care must be taken when adopt- 
ing quoted results from previous studies. For example, the deter- 
mination of the distance to water masers is based on the assump- 
tion of Keplerian motions. Furthermore the pulsation of Cepheids 
relies on the star's self-gravitation, therefore unusual gravitational 
behaviour could lead to a modification of their period-luminosity 
relation. However there are currently no signs of anomalous be- 
haviour in the data I jJain et al.|2012| l. 



3.2 Baryon Acoustic Oscillations 

Encoded within the distribution of galaxies is a useful comoving 
ruler, a modest excess of galaxy pairs at ~150 Mpc, relating to 
the distance sound waves were able to propagate prior to recom- 
bination. In Fourier space this feature is referred to as the Baryon 
Acoustic Oscillations (hereafter BAO). For a given BAO observa- 
tion, the relevant cosmological information can often be distilled to 
a single datapoint, that of Dy{z)/rs, where rs is the sound hori- 
zon at the baryon drag epoch, and the volume element Dv(z) is 
defined in terms of the angular diameter distance Da{z) and Hub- 
ble parameter H{z) 



Dv{z) = 



{l + zfDi{z) 



H{z 



(10) 



The most conservative treatment of BAO data would be to treat 
Ts as an unknown quantity, and eliminate it by taking the ratio of 
Dv{z) at two separate redshifts. However since the CMB provides 
us with a good understanding of the conditions in the early Uni- 
verse, we have the opportunity to determine Vs for a given set of 
cosmological parameters. For our purposes this quantity is derived 
using equation (6) from |Eisenstein & Hu| ( |1998] l. 

Where stated, we utilise measurements of Dy{z)/rs at two 
different redshifts. [Anderson et al.|p012[ > present BAO constraints 
from the Baryon Oscillation Spectroscopic Survey (BOSS), and 
we also include a measurement from the Sloan Digital Sky Survey 
(SDSS) Luminous Red Galaxies (Padm anabhan et al.|201 2). Both 
studies present results before and after a reconstruction method has 
been applied, a technique which increases the prominence of the 
BAO by partially reversing the effects of gravitational infall. We 
conservatively adopt the BAO results prior to reconstruction, since 
we wish to avoid making assumptions relating to the growth of 
density perturbations. The change in our results is negligible, how- 
ever, when the reconstructed values are included. For the BOSS 



CMASS sample, the 'consensus' result ^Anderson et al.|2012 l) aver- 
ages the power spectrum and correlation fiinction analyses, giving 
Dv/niz = 0.57) = 13.67 ± 0.22. The SDSS datapoint is given 
by D\r /rs{z = 0.35) = 8.89 ± 0.31. We neglect any covariance 
in the two data points generated by the small overlap between the 
SDSS and BOSS volumes, which spans less than 9% of the BOSS 
volume ( An derson et al.|20l2j . 

It should be noted that when combined with WMAP7, recent 
BAO results suggest values of Ho which are in mild tension with 
the rneasure ment from |Riess et al.| ( |2011[ l. For example [Anderson] 
8.4 ± 1.3kms"^ Mpc"\ while stud- 



et al. 



2012) find Ho 



ies from 6dFGS and SDSS measure 67± 3.2kms~^Mpc~^ and 



69.8 ± 1 .2 km s ^ Mpc ^ respectively i Beutler et al. 



2011 



Mehta 



[et al.|201^ . We present separate results generated using the Ho and 
BAO data in Section|6] 



3.3 Cosmic Microwave Background {£ > 100) 

We utilise the temperature Cf^ and polarisation Cf^ spectra from 
the WMAP 7-year data ^Larson et al.|201 T \. This serves two pur- 
poses, first it offers strong constraints on a number of cosmological 
parameters, but it also serves to directly assist in the measurement 
of the growth of structure. Our current position in space and time, 
z — 0, acts as the benchmark from which all cosmological dis- 
tance measurements are made. However when mapping the growth 
of structure, it is the epoch of recombination at z ~ 1100 which 
effectively acts as our starting point. The CMB provides a precise 
measure of the amplitude of density fluctuations in the early Uni- 
verse. Combining this with a measure of the clustering amplitude 
at low redshift provides a strong lever for constraining the growth 
rate between these two epochs. 

Large angular scales of the CMB are sensitive to modified 
gravity models, and these are treated separately in j ]4.3[ The im- 
plications for modified gravity have already been studied (see for 
example fZuntz et al.|2012^ , and here our primary focus is to explore 
the constraints available from lensing. For the range of values con- 
sidered in this study, and Eq were found to impose negligible 
change in either Cf^ or Cf^ beyond the £ = 100 multipole (see 
Section [43| >. Unless otherwise stated we simply truncate the largest 
angular scales, such that only £ ^ 100 are included in our analysis. 
In doing so we negate the impact of the ISW, and utilise WMAP7 
only to assist in measuring the conventional cosmological parame- 
ters. For further details of our treatment of the CMB, see Appendix 

El 



4 GRAVITATIONAL DATA 
4.1 Cosmic Shear 

Gravitational lensing by large scale structures in the Universe co- 
herently distort the observed shapes of distant galaxies. This cos- 
mic shear distortion can be used to directly probe dark matter, and, 
when analysed in tomographic redshift bins, can also reveal the 
growth of structure. It also acts as a useful tool for testing the laws 
of gravity (White & Kochanek 200T][Dore et al.|2007^ . We use to- 
mographic cosmic shear data from CFHTLenS, published in [Ben-[ 
[jamin et ar] ( [2012} . CFHTLenS is a 154 square degree deep multi- 
colour survey optimised for weak lensing analyses, where the data 
were obtained as part of the CFHT Legacy Survey on the 3.6m 
CFHT. Weak gravitational lensing is widely recognised as one of 
the most powerful probes of cosmology, but also one of the most 
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Figure 2. The cosmic sliear data from CFHTLenS. Eacli panel contains the two lensing correlation functions, £,+ {d) (filled) and ^—(9) (empty) with the 
C-(^) points slightly offset for clarity. The source galaxies are divided into two redshift bins, a 'low' redshift bin with 0.5 < z ^ 0.85 and a 'high' 
redshift bin with 0.85 < 2 ^ 1.3. The left, middle and right hand panels display data from the low-low, low-high, and high-high correlations. In each case, 
the black lines are the set of theoretical curves corresponding to 5+ (9) (upper) and 5_ (S) (lower) for the best-fitting ACDM cosmology. The dashed lines 
correspond to a cosmological model with the same expansion histoiy H(z) as the best-fitting model, but with a boost in the effective gravitational constant 
(/^O = 0.5; So = 0). The dot-dashed lines demonstrate the signal one expects if the present-day (z = 0) gravitational potentials are associated with Newtonian 
gravity (/.jq = 0; Sq = -0.5). 



technically challenging, as the weak cosmic shear distortion must 
be separated from other sources of image distortion induced by the 
atmosphere, telescope and detector. The CFHTLenS analysis rep- 
resents the current state-of-the-art in weak lensing data processing 
( |Erben et al.|2009| l, weak lensing measurement (Miller et al.|2012| l, 
photometric redshift measurement fHildebrandt et al. 2012 1, and 
systematic error analysis ( [Heymans et al. 2012) . Benjamin et al.| 
( |2012^ also use an angular cross-correlation analysis to verify the 
accuracy of the measured redshift distributions that we use in this 
analysis. These results, in combination with the cosmology insensi- 
tive joint shear-redshift systematic test presented in |Heymans et al.| 
( |2012) , allows us to conclude that the analysis of the survey data is 
robust within our uncertainties and residual systematics levels and 
can therefore be used to test the laws of gravity. 

[Benjamin et al.| ^2012^ analyse two high-redshift bins with 
0.5 < z ^ 0.85 and 0.85 < ^ ^ 1.3, presenting the auto- 
correlation and cross-correlations of the angular two-point shear 
correlation functions Cl{0) over scales 1 < ^ < 40 arcmin, as 
shown in Figure |2] The choice to use only two broad high red- 
shift bins in this tomographic analysis was motivated by the de- 
sire to mitigate the impact of intrinsic galaxy alignments on the 
tomographic measurements (see for example 'Heym ans et al.|2006l 
|Joachimi et al .|20TT. and references therein). For this combination 
of high redshift bins, the contamination to the signal is expected to 
be at the percent level and so we follow |Benjamin et al.| ( |2012) l by 
neglecting this source of error in our analysis. We use a covariance 
matrix between the resulting correlated data points as estimated 



from N-body lensing simulations ( Hamois-Deraps et al.|2012 1, dis- 
cussed and verified in |Kilbinger et al. i2012i and Benjamin et al.| 
([2012). 

Theoretical predictions for these correlation functions are gen- 
erated using the relation 



j_ r 
2^/0 



(11) 



where P^'-' {£) is given in equation jsj and J0/4 are the zeroth and 
fourth order Bessel functions. 

On the smallest angular scales studied by cosmic shear it is 
likely that baryonic physics will disturb the density perturbations 
( [Semboloni et al.|2011^ . This is also the regime where our prescrip- 
tion for modelling the non-linear correction to the matter power 
spectrum (see below) begins to break down. Due to the very dif- 
ferent nature of the zeroth and fourth Bessel functions, the ^_ (6) 
correlation function is sensitive to considerably smaller scales than 
^+(0) for a given 9 value. For the ^_ statistic we therefore only 
include the larger angular scales 9 > 10'. These data points are 
presented in Figure |2] alongside theoretical predictions from the 
best-fit ACDM cosmology and two illustrative modified gravity 
scenarios. From left to right, the three panels in Figure |2] relate 
to the correlations of the low-low, low-high, and high-high redshift 
bins. A more detailed justification for the angular cuts and their 
implications are discussed in [Benjamin et al.1 ( |20l2^ ; in any case 
we find no significant change in our results when the small angular 
scales of ^_ {9) are included due to the high correlation between 
the points and their relatively low signal-to-noise. We further in- 
vestigate the impact of truncating smaller length scales in Kitching 
et al (in prep). 



4.2 Redshift Space Distortions 

Doppler distortions of the galaxy clustering pattern provide an in- 
dependent probe of gravitational physics, by revealing the coher- 
ent flows of galaxies which trace the growth of cosmic structure. 
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The principal effect on large scales is an amplification of appar- 
ent galaxy clustering along the line-of-sight, breaking the statistical 
isotropy. This anisotropy can be extracted from the 2-point cluster- 
ing statistics of a galaxy redshift survey, yielding a measurement 
of f{z)as{z). Here / is the growth rate of structure, expressible 
in terms of the growth factor D{a) as / = ciln_D/cilna where the 
growth factor describes the evolution of the amplitude of a single 
linear density perturbation 5{a) = D{a)5{l). The function as{z) 
quantifies the amplitude of the linear clustering of matter at a given 
redshift, and is defined as the r.m.s. density variation in spheres of 
comoving radius 8 h~^Mpc. The degeneracy between / and erg 
arises due to the unknown galaxy bias factor 6: the observables are 
the anisotropy in the galaxy clustering pattern which depends on 
j3 = f /b, and the galaxy clustering amplitude which depends on 
basiz). 

Measurements of redshift-space distortions from a galaxy red- 
shift survey are degenerate with a second isotropy-breaking mecha- 
nism: the Alcock-Paczynski distortion ( |Alcock & Paczynski|1979[ l, 
which arises if the fiducial cosmological model we adopt for con- 
verting redshifts and angles to physical distances differs from the 
true cosmological model (Ballinger et al. 1996; Matsubar a & Suto| 
[T996)|Seo & Eisenstein|2003,. Simpson & Peacock.2010j . The am- 
plitude of Alcock-Paczynski anisotropy depends on the quantity 




Figure 3. The growth data from WiggleZ jBlake et aI.|20T2 ll and 6dFGS 
(Beutler et al . 20121. The black line is the corresponding theoretical curve 
for the best-fit ACDM cosmology. The dashed lines correspond to the same 
cosmology but with a boost in effective gravitational constant fiQ . 



F{z) = {1 + z)Da{z)H{z)/c, 



(12) 



in terms of the angular diameter distance Da{z) and Hubble ex- 
pansion parameter H{z). Although measurements of F and /erg 
from a galaxy redshift survey are strongly correlated, the angular 
distortions imprinted in the clustering pattern by these two effects 
are sufficiently different that the degeneracy can be broken. 

The measurements of (/erg, F) that we use in this study are 
derived from the WiggleZ Dark Energy Survey published in |Blake| 
|et al. ^20 12 1 and from 6dF Galaxy Survey (6dFGS), as presented 



Beutler et al.|p012l l. Figure [3] compares the measurements of 



fas from WiggleZ (squares) and the 6dFGS (diamond) to the best- 
fit ACDM cosmology and two example modified gravity scenarios 
with fio = 0.5 and /io = 1.0 (dashed). 

The WiggleZ Survey ( |Drinkwater et al.|[20T0l > at the 3.9- 
m Anglo- Australian Telescope has mapped a cosmic volume ~1 
Gpc"^ over the redshift range < z < 1. By covering a total of 
~800 deg'^ of sky the WiggleZ survey mapped about 100 times 
more effective cosmic volume in the z > 0.5 Universe than previ- 
ous galaxy redshift surveys. Target galaxies were chosen by a joint 
selection in UV and optical wavebands, using observations by the 
Galaxy Evolution Explorer satellite (GALEX) ( jMartin et al.|2005^ 
matched with ground-based optical imaging from the Sloan Digital 
Sky Survey ( [York et al.|2000j l in the Northern Galactic Cap, and 
from the Red Sequence Cluster Survey 2 (RCS2) jGilbank et al.| 
|201 l| l in the Southern Galactic Cap. A series of magnitude and 
colour cuts ( [Drinkwater et al.|2010| l were used to preferentially se- 
lect high-redshift star-forming galaxies with bright emission lines, 
that were then observed using the AAOmega multiobject spectro- 
graph ( [Sharp et al.|2006^ in 1-hr exposures. 

The 6dFGS is a near-infrared selected combined redshift and 
peculiar velocity survey covering nearly the entire southern sky 
with Galactic latitude |&| < 10°. It was undertaken with the Six- 
Degree Field (6dF) multi-fibre instrument on the UK Schmidt Tele- 
scope from 2001 to 2006 ( jjones etal .2Q04 20091. The redshift sur- 
vey covers 125,000 galaxies with 8.75 ^ ^ 12.75 and median 
redshift z = 0.052. 

In this study we use the WiggleZ joint measurements of 



(/(Tg, F) in three different redshift slices with effective redshifts 
z = 0.44, 0.60, and 0.73. The marginalised measurements of 
the growth rate in these redshift slices are fas{z) = (0.41 ± 
0.08, 0.39 ± 0.06, 0.44 ± 0.07). The marginahsed measurements 
of the Alcock-Paczynski distortion parameter are F = (0.48 ± 
0.05,0.65 ± 0.05,0.86 ± 0.07). The correlation coefficients be- 
tween these pairs of values are (0.73, 0.74, 0.85). These measure- 
ments were derived by fitting the non-linear clustering model pro- 
vided by |Jennings et al.| ( |201 l| l to the 2D WiggleZ power spectrum 
over the range of wavenumbers < k < 0.2 h Mpc^^, marginal- 
ising over the linear galaxy bias factor and a small-scale velocity 
dispersion (Ty. Blake et al.| ( (2011[ l carried out a suite of systematics 
tests to demonstrate that the derived results did not depend on either 
the fitting range or the non-linear model of redshift-space distor- 
tions adopted; we refer the reader to this paper for further details. 
There is some overlap (~500 sq. degrees) between the sky areas 
surveyed by WiggleZ and the SDSS/BOSS surveys that are used 
in our results which are stated to include BAO. However, based 
on the nature of the geometric BAO data being very different to 
that used for the RSD gravitational measurements we assume that 
any covariance in these derived results is negligible. In essence, 
when considering the redshift-space galaxy power spectrum, the 
location of the acoustic peak bears little connection to the power of 
the quadrupole. 

We also include the single datapoint ,f(Tg,[z = 0.067) = 
0.423 ± 0.055 ( Beutler et al. 2012 ) from the 6dFGS. Due to its 
low redshift, this result benefits from having negligible sensitivity 
to the Alcock-Paczynski distortion, and a large sensitivity to our 
modified gravity parameter /io, as can be seen in Figure |3] The 
6dFGS measurement was derived using the non-linear model of 
Jennings et al. (2011), including wide-angle corrections, fit to the 
2D galaxy redshift-space correlation function for transverse separa- 
tions rp > 16 Mpc. The growth-rate measurements obtained, 
marginalized over bias, were consistent with fitting an empirical 
streaming model to the results (see |Beutler et al. '2012 for further 
details). While the WiggleZ and 6dFGS analyses differ in their util- 
isation of the power spectrum and correlation function respectively, 
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Figure 4. Modifications to the temperature (upper panel) and temperature- 
polarisation (lower panel) power spectra of the CMB. In each panel the 
black line is the corresponding theory curves for the best-fitting ACDM 
cosmology. The dashed and dash-dot lines represent the same modifications 
to our modified gravity parameters as seen in Figure|2] Superposed on these 
theoretical predictions are the binned data points from the 7— year WMAP 
observations jLarson et al.|201 f) . 



their approach to modelling non-linear corrections and galaxy bias 
remains consistent. 



4.3 Cosmic Microwave Background {I < 100) 

When in the presence of a time-dependent gravitational potential, 
a particle may extract or lose energy from its surroundings. It is 
this mechanism which space probes exploit when passing close to 
planets for a gravitational slingshot. For CMB photons, the gravita- 
tional source is associated with large scale density perturbations. In 
GR, if the photon passes through a high density region of the Uni- 
verse, the dark energy-induced decay of the gravitational potentials 
$(a;, t) and '^{x,t) results in a net energy gain, with a blue-shifting 
effect. Conversely, passing through a void results in a net loss of en- 
ergy. The resulting temperature perturbation may be expressed by 
the line integral 



where r) is the conformal time defined by Ar) = At /a, and the inte- 
gral follows the path of the photon. We can see from equation l |13[ > 
that this ISW effect is another probe which, like gravitational lens- 
ing, retains sensitivity to both gravitational potentials. However the 
ISW signal within the temperature anisotropies is only expected to 
be significant on large angular scales, I < 100 as illustrated by 
Figure|4]which compares the temperature (upper) and temperature- 
polarisation (lower panel) GR power spectra of the CMB (solid) 
with data points used in this analysis, from the 7— year WMAP 
observations ( [Larson et al.||20lT| . The dashed and dash-dot lines 
represent the same modifications to our modified gravity param- 
eters as seen in Figure |2] with fio ~ 0.5; So — (dashed) and 
Ho = 0; So — —0.5 (dot-dashed). These modified gravity models 
of the CMB power spectrum deviate from GR only on the largest 
angular scales. Observational ISW studies are therefore heavily re- 
stricted by cosmic variance, and future experiments cannot be ex- 
pected to improve upon WMAP, unless we begin to include cross- 
correlations with large-scale structure. Note that throughout this 
work, we use the term ISW to refer to the phenomenon which arises 
in the CMB power spectrum, not the cross-correlation between the 
CMB and the distribution of galaxies. 

For any metric theory of gravity, density fluctuations on scales 
beyond the horizon are subject to a stringent consistency relation. 
Since material cannot be transported at superluminal velocities, su- 
perhorizon perturbations are destined to evolve along a path deter- 
mined by the global expansion ( [Berts chinger 2006). Therefore the 
constraints presented which incorporate the ISW should be treated 
with this in mind, as the ISW involves scales very much larger than 
those studied by gravitational lensing, a regime where there is a 
stronger expectation that the values of /lo and Eo are close to zero. 



5 METHODS 

Using a combination of geometric data (presented in f|3j and gravi- 
tational data (presented in !|4j, we use a maximum likelihood analy- 
sis to constrain our phenomenological modified gravity parameters 
^0 and So. 

To generate the likelihood values we modified CosmoPMC 
( [Kilbinger et al.|[20IH , an adaptive importance sampling code, 
to incorporate the alterations to the growth of structure (equa- 
tion |7j and cosmic shear signal (equation [5}. A custo mised 
version of MGCAMB jHojjati et ar]|20TT| [Lewis et al.1[2000l > 
was used to evaluate the CMB anisotropies for a given set 
of cosmological parameters, which are then assessed with the 
WMAP likelihood cod^ Our analysis spans the parameter 
set {Qm, i^b, h,T,nB, As z , ^o,'^o}, where r is the optical 
depth. Us the spectral index, Asz the amplitude of the Sunyaev- 
Zel'dovich template for the CMB ( [Komatsu & Seljak|2002l[Spergel[ 
et al.|20()7 1, and controls the primordial amplitude of matter 
perturbations. Where specified we also marginalise over a time- 
independent effective equation of state for dark energy w and 
global curvature SIk, in addition to the aforementioned set of pa- 
rameters. We adopt flat priors throughout, spanning the following 
parameter ranges: e [0.05; 1], fib G [0.01; 0.1], h £ [0.5; 1], 



T e [0.02; 0.3], G [0.8; 1.2], Asz G [0; 2], G [1.8; 3.5] x 



^ http://lambda.gsfc.nasa.gov 



© 0000 RAS, MNRAS 000, 000-000 



8 F. Simpson et al. 



10-^ MO G [-4.0; 4.0], Eq G [-4.0; 4.0], w G [-1.5; -0.3], and 
f^K G [—0.1; 0.1]. In most cases these boundaries extend well be- 
yond the regions of high likelihood, except for Asz which is not 
well constrained, however its impact on the other parameters is neg- 
ligible. In order to generate visualisations of the generated samples 
we make use of code adapted from CosmoloGUP] 

In generating the matter power spectrum, the CosmoPMC 
package makes use of the [Eisenstein & Hu| ( |1998[ l transfer func- 
tion. The growth factor is determined in accordance with the mod- 
ified equation JTl, before adding a nonlinear correction using the 
prescription of [Smith et al.[(2003] l. To compensate for cosmologies 
where w 7^ —1, a further modulation of the nonlinear power is 
required, as presented in [McDonald et alT]i 2006 1. We make use of 
the scaling correction provided in Schrabback et al.[ ( |2010^ to apply 
this formalism across our full range of parameter space. However 
it should be noted that our results only have a weak dependence on 
this correction. The mean values of /.to, Eq, and w all change by 
less than 0.02 when the prescription from |Schrabback et al.[j2010) 
is omitted, which is well within our statistical errors. The value of 
the Eo parameter has no impact on structure formation, and while 
we may expect nonzero values of /xq to alter the nonlinear correc- 
tion to the matter power spectrum, we do not attempt to correct 
for it here. We expect the prescription of [Smith et al.[ ( [2003| > to re- 
main a reasonably good approximation in the presence of a simple 
scale-independent modification to the growth of structure such as 
that induced by /iq. To first order the correction associated with 
nonlinear power is simply a function of the ambient linear power 
at that epoch; how quickly it arrives at this state is of less impor- 
tance. In any case, modified behaviour on nonlinear scales would 
still weaken the fit to the data achievable within the framework of 
GR. 



6 RESULTS 
6.1 ACDM 

We begin with the simple case of a flat ACDM geometry. In this 
scenario we are making a strong assumption about the form of 
the cosmic expansion history H{z), but in doing so make signif- 
icant gains in the precision of our likelihood contours. Any appar- 
ent deviation from GR witnessed at this stage must be interpreted 
with caution due to the strong assumptions of this relatively simple 
model. 

The left panel of Figure [3] illustrates the variety of 1- and 2~ 
a confidence intervals for different combinations of datasets. Each 
of these contours include the auxiliary data of WMAP7 (£ ^ 100) 
and a prior on Hq. The horizontal (green) contours reflect the pure 
Ho measurement established with the WiggleZ data. The Eo pa- 
rameter is unconstrained since the peculiar motions of galaxies are 
slow (n <C c), and are therefore insensitive to the spatial curvature 
<l>(a;, t). The near-vertical (red) contours correspond to the cosmic 
shear data from CFHTLenS. Unsurprisingly, this is predominantly 
a measure of Eo due to its direct influence on the lensing poten- 
tial. The mild sensitivity to fio arises from the amplitude of the 
matter power spectrum entering in equation A positive value 
of fj,o enhances the growth of structure and thus strengthens the 
cosmic shear signal (for a given primordial amplitude). This can 
be offset with a negative value of Eo, weakening the deflection 
of light and leading to the negatively-sloped degeneracy direction 



http://www.sarahbridle.net/cosmologui 



seen in the red contours of Figure [5] When combined (blue), the 
near-orthogonal degeneracy directions of WiggleZ and CFHTLenS 
produce a considerable reduction in the permitted area of param- 
eter space. The tightest constraints (white) emerge from adding in 
BAO from BOSS and the ISW large scale anisotropics {£ < 100) 
from WMAP7. None of these combinations show any preference 
for non-zero values of either /io or Eo, suggesting that GR remains 
a successful model of gravitation on cosmological scales, for both 
relativistic and non-relativistic particles. 

The mean and l-a errors in the modified gravity parameters 
are found to be fio = 0.05 ± 0.25 and Eo = 0.00 ± 0.14. The 
results for other combinations of datasets can be found in Table [Tj 

How sensitive are we to the choice of auxiliary data? The right 
hand panel of Figure[5]contains the same datasets and layout as the 
left hand panel, except here we have replaced the prior on Hq with 
the BAO data summarised in Section [T2| The BAO prefer a higher 
value of Q,rn (or equivalently, a lower value of Hq), and this does 
push the lensing data to predict a lower value of erg, which pushes 
the bulk of the likelihood towards negative values for Eo. This re- 
sults in constraints of /io = -0.12±0.23andEo = -0.17±0.10. 
The combination of CFHTLenS and RSD data (blue) can be com- 
pared with the constraints derived from the full WMAP7 power 
spectra, including the ISW, and BAO (yellow). 



6.2 uiCDM 

The dangers of confusing the signatures of modified geometry and 
modified gravity have previously been highlighted in the literature 
(see for example [Simpson & Peacock|2010[ l. In order to claim that 
we are testing GR and not simply ACDM we must therefore al- 
low greater freedom in the cosmic expansion history H{z). This is 
inevitably accompanied by some loss of precision, so we shall re- 
strict ourselves to adding one extra parameter. In Figure|6]we utilise 
both Hq and BAO, in combination with CFHTLenS, RSD and the 
small angular scales of WMAP7, conservatively omitting the ISW 
information for the reasons outlined in Section [43] While the green 
contour in the left panel corresponds to the standard flat ACDM 
background, the blue contours presented in the left and centre pan- 
els allow the effective dark energy equation of state w to deviate 
from —1. The centre panel illustrates a degeneracy which causes 
the confidence contours to broaden. More positive values of w lead 
to a greater prevalence of dark energy at higher redshift (for a given 
value of flm) and thus a slowing of the growth of structure. This can 
be compensated by enhancing the value of fi, which accelerates the 
growth of structure, leading to a similar growth history f{z) to GR. 
For a wCDM cosmology the mean and 68 per cent confidence in- 
tervals are given by /xo = -0.59 ± 0.34, Eo = -0.19 ±0.11, and 
w — —1.19 ± 0.1. More detailed results are presented in Table[T] 



6.3 owCDM 

We further generalise the expansion history by allowing for non- 
zero values of the global curvature flk, using the same combination 
of data sets as in j j6.2| The red contours in all three panels of Figure 
[6] demonstrate the relatively small impact that non-flat geometries 
have on our confidence limits. The mean and 68 per cent confidence 
intervals are given by /io = -0.65 ± 0.34, Eo = -0.26 ± 0.12, 
Qk = 0.096 ± 0.087, and w = -1.16 ± 0.10. 
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Figure 5. Left: Constraints on the modified gravity parameters in a flat ACDM background from redshift space distortions (green), weak lensing (red), and 
combined (blue). (68 and 95 per cent CL). The dashed and solid contours represent the 68 and 95 per cent condence intervals respectively. Two auxiliary 
datasets are used here to break degeneracies with the conventional cosmological parameters. These are the small-angle anisotropics from WMAP7 (£ ^ 100), 
and a prior on Ho from jRiess et al.|20TT) . The cross positioned at the origin denotes the prediction of General Relativity. Right: The red, green and blue 
contours are the same as the left panel, except the prior on Ho has been replaced by measurements of the Baryon Acoustic Oscillations as detailed in Section 
|3.2| The yellow contours signify the constraints derived from the full WMAP7 power spectra, including the large angular scales {£ < 100). The white contours 
in the left hand panel show the constraints when all data sets are analysed in combination. 



Table 1. Parameter constraints for different combinations of dataset and parameter space. The three different backgrounds we explore are flat ACDM, flat 
uiCDM and the non-flat uiCDM which we denote as oiiiCDM. All constraints make use of the small scale anisotropics from WMAP7 {£ ^ 100), while those 
marked ISW also utilise the larger angular scales. 



Background 


Dala CMB + Hq ) 




So 






^8 






ACDM 


CFHTLenS 


L2 ± 2.2 


-0.17 ± 0.28 


0.260 ± 0.024 


71.3 ± 2.3 


0.96 ± 0.24 






ACDM 


RSD 


0.08 ± 0.25 




0.254 ± 0.023 


71.8 ± 2.2 


0.801 ± 0.024 






ACDM 


CFHTLcLiS + RSD 


0.05 ± 0.25 


0.00 ±0.14 


0.256 ± 0.023 


71.6 ± 2.2 


0.804 ± 0.023 






wCDM 


CFHTLenS + RSD + BAO 


-0.59 ± 0.34 


-0.19 ± 0.11 


0.272 ± 0.015 


71.7 ± 1.8 


0.820 ± 0.025 


-1.19 ± 0.10 




oxL-CDM 


CFHTLenS + RSD -I- BAO 


-0.65 ± 0.34 


-0.26 ±0.12 


0.289 ± 0.021 


71.6 ± 1.8 


0.833 ± 0.028 


-1.16 ± 0.10 


0.0096 ± 0.0087 


Background 


Data f-i- CMB + BAO) 
















ACDM 


CFHTLenS 


].0 ± 2.2 


-0.31 ± 0.24 


0.305 ± 0.018 


67.2 ± 1.4 


0.95 ± 0.24 






ACDM 


RSD 


-0.13 ± 0.23 




0.300 ± 0.017 


67.5 ± 1.4 


0.807 ± 0.023 






ACDM 


ISW 


-0.62 ± 0.64 


0.34 ± 0.22 


0.294 ± 0.015 


68.2 ± 1.2 


0.768 ± 0.068 






ACDM 


CFHTLenS + RSD 


-0.12 ± 0.23 


-0.17 ± 0.10 


0.300 ± 0.017 


67.6 ± 1.4 


0.807 ± 0.023 






ACDM 


CFHTLenS + RSD + ISW + Hq 


-0.06 ± 0.21 


-0.010 ± 0.068 


0.273 ±0.010 


70.10 ± 0.97 


0.803 ± 0.022 







1 ALTERNATIVE PARAMETERISATIONS 7.1 7 Parameter 

A widely adopted parameterisation for quantifying anomalous 
structure growth is the growth index 7 ( [Wang & Steinhardt|1998[ 
|Linderp005l > 

In this section we consider the relationship between the E) pa- 
rameterisation adopted in tiiis work, and thtree other popular param- 

eterisations. /(z) ~ , (14) 
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Figure 6. In this figure we extend our analysis to more general expansion histories. All of the above contours use a combination of cosmic shear and redshift 
space distortions, combined with the geometric constraints from Hq, BAO, and WMAP7 {I ^ 100). The green set of contours in the left panel corresponds 
to the background expansion of flat ACDM. For the blue contours in the left and centre panels we allow the effective dark energy equation of state to deviate 
from to = — 1, highlighting a significant degeneracy in the w — no plane. Finally, the red contours appearing in each of the panels allow non-zero curvature 
in addition to the toCDM background. 




Figure 7. Confidence limits in the 7 - Qrn plane, as derived from red- 
shift space distortions (WiggleZ + 6dFGS) in green, and cosmic shear 
(CFHTLenS) in red (68 and 95 per cent CL), using small-angle anisotropics 
from WMAP7 (£ ^ 100), and a prior on Hq as auxiliary data sets. Note that 
here 7 is the only modified gravity parameter, as we have set S(a) = 0. 



where setting 7 = 0.55 offers an estimate of the growth factors 
in wCDM cosmologies with sub-percent accuracy. Equation (23) 
from |Linder & Cahn| ( (2007l l gives us the relation 



3 1 



7 



M(a) 
l-n„,(a) 



5 — 6ui 



(15) 



for a constant equation of state w. In a flat ACDM Universe, this 
parameter may be expressed as a simple function of fio and f^A, 



7 = 



6^ 

n 



JM)_ 



(16) 



Applying the prescription for the growth rate in equation ^14) 
with a constant value the 7 parameter does not offer a perfect re- 



Figure 8. The same combinations of data as the left hand panel of Figure 
|5] except we have replaced the /jq parameter with 7 using the relation given 
by equation \16) . 



production of the growth of density perturbations in GR. This rela- 
tionship provides a useful means to derive constraints in terms of 7 
without compromising the precision of the calculation. Constraints 
on this quantity are presented in the 7 — Qm plane in Figure |7] 
Gravitational lensing is not uniquely defined when using the 7 pa- 
rameter alone, so here we have simply set So = such that the 
lensing potential is unaltered from its conventional form. We find a 
value of 7 = 0.52 ± 0.09 for the case of a ACDM background. In 
Figure[8] we then let Eo vary, presenting constraints in the 7 — Eo 
plane. 



7.2 Eg Parameter 



Proposed by [Zhang et al.|p007[ l, the Eq parameter encompasses 
contributions from both gravitational lensing and the motions of 
galaxies. This is in the form of a ratio which by construction nulli- 
fies contributions from both linear galaxy bias and the amplitude of 
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Figure 9. Confidence Hmits on the So and /jq parameters derived from the 
measurement Eq{z = 0.32) = 0.392 ± 0.065 ^Reyes et al.|2010| . The 
inner blue contours are the same as those in Figure [5] Here we assume a 
ACDM cosmology. 




Figure 10. The same combinations of data as the left hand panel of Figure 
|5] except here we have replaced the Eg parameter with the gravitational slip 
(^0 as defined in equation jl9) . 



the matter power spectrum. The result is a quantity which is sensi- 
tive to both the lensing potential and the rate of structure forination 



Eg 



3//2(l + z)/35 



(17) 



where jS = f/b. This technique has been applied to SDSS data 
in [Reyes et al.| ( (2010^ . We cannot directly measure Eg with the 
CFHTLenS and WiggleZ datasets, since they do not overlap in sky 
coverage. However we can perform a direct mapping of the result 



from Reyes et al. (2010i onto the {fi, E) plane, utilising the rela- 
tionship 



Er 



(18) 



where the dependence on fi arises through its influence on f{z), 
the linear growth rate at the mean redshift of the survey. The results 
are presented in Figure[9] where the red contours correspond to the 
confidence limits, in the fio — Eq plane, arising from a combination 
of the measurement Eg {z = 0.32) = 0.392 ± 0.065 ([Reyes et al 



[2010| > with the BAO data from Section [T2| which provide a con- 
straint on Qm. The CMB does not provide further assistance in this 
case because the Eg parameter is insensitive to the amplitude of 
perturbations. For comparison, the inner blue contours match those 
in the right hand panel of Figure [5] We note that GR can readily 
generate departures from the conventional Eg value. For exam- 
ple in equation the ratio Qrn/f{z) changes when S7fe 7^ or 
U) 7^ — 1. Similarly, radically different modified gravity models can 
generate an Eg value which appears consistent with the ACDM 
value. This is apparent from the degeneracy contour seen in Figure 
|9] where large values of both jj, and E enhance the growth rate and 
lensing signal, leaving their ratio in equation l [17[ ( unchanged. 



7.3 Gravitational Slip 

Rather than modulating the sum of the two potentials using S(a), 
it is often convenient to work in terms of the ratio of the potentials. 
One definition in this form is known as the gravitational slip 



C = 1- 



(19) 



which we expect to be zero in GR, in the absence of anisotropic 
stress. This has been the focus of much theoretical exploration, for 
example [Ferreira & Skordi^ j2010| l demonstrate that for second- 
order theories we generally expect to find this value scales with the 
Newtonian potential and its time derivative. 

Recasting in terms of (/j., E) using equations ^ and |5]l 



C(«) = 



2[E(a)-M(a)] 



(20) 



1 + 2E(a) - ^(a) ' 

As an illustration in Figure [lO] we reformulate the x-axis from Fig- 
ure[5]in terms of (^0, the gravitational slip evaluated at 2 = 0. One 
striking feature of these contours are the long tails out to nega- 
tive values of and fi, corresponding to where the denominator 
<1> ventures close to zero. This suggests that the slip parameter may 
possess a pathological likelihood surface, and so for the purposes of 
performing Monte Carlo evaluations these contours are more safely 
generated as a derived parameter. 



7.4 Gravitational Potentials 

The redshift-space distortion data, used in this analysis, exhibits 
a peak sensitivity to a modified gravity scenario at approximately 
~ 0.5 (see Section [8|. In Figure [TT] we therefore consider frac- 
tional deviations in the two gravitational potentials at this redshift, 
relative to their GR prediction for the same mass distribution. 
These are given by 



A* 



(21) 



where AvE" = — $gr. For our conservative combined lensing, 
RSD, Ho and high i WMAP7 data analysis we find cosmological 
deviations in the Newtonian potential and curvature potential from 
the prediction of General Relativity to be A*/* = 0.05 ± 0.25 
and A$/$ = -0.05 ± 0.3 respectively (68 per cent CL). 
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Figure 11. Here we explore fractional deviations in the two gravitational 
potentials, the Newtonian potential ^ and the curvature potential <1?, from 
the GR value at z = 0.5. The prescription for this is given by equation 
The contours represent the same combinations of data as those in the left 
hand panel of Figure[5] 



8 THEORETICAL MODELS 

Below we briefly review some of the theoretical models which 
could generate a departure from /io = So = 0, and interpret the 
implications of our results. There is such a plethora of modified 
gravity models, that no single choice of parameterisation can ade- 
quately encompass all of them. This is a situation reminiscent of the 
dark energy equation of state, w{z), except here we are faced with 
uncertainty not only in the functional form of the time-dependence, 
but also in its scale-dependence. So how can we relate a given 
(/X, E) constraint to a specific model? The observed parameters /io 
and Eo may be interpreted as a weighted integral over the true func- 
tional form jj,{k,z), such that 



(t){k, z) lJ.{k, 



Qa{z] 



dkdz . 



(22) 



If we perform a scale- and time-dependent principal component 
analysis (see for example Zha o et al.|2009 '), then the weight func- 
tion <j){k, z) may be expressed in terms of the principal components 
ei{k, z) and the errors associated with their corresponding eigen- 



values o-{aj) I Simpson & Bridle 2006 1, such that 



(t){k,z) = 



ei{k, z) JJ e^jk', z')dk' dz'/a^{ai) 
E, [JJe,ik",z")dk"dz'f /aHa,) 



(23) 



The analysis of redshift space distortions in |Blake et al.| ( [20T2) in- 
cludes information from the galaxy power spectrum up to a max- 
imum wavenumber fcmax ~ 0.2^Mpc~^, corresponding to the 
regime over which the density and velocity fields are sufficiently 
linear for our theoretical models to remain valid. Since the number 
of Fourier modes increases towards higher k, the scale-dependent 
component of 0(fc, z) peaks close to this value of fcmax, and 
(f){k, z) = for k > fcmax. We evaluate the redshift-dependence 
of the weight function (j){z) associated with the combined WiggleZ 
and 6dFGS data of Figure [3] following the prescription of |Simp-| 
|son & Bridle| l|2006), and this is shown to peak at z ~ 0.5 as il- 
lustrated in Figure [T2] In the following subsections we utilise the 
weight function 4>{z) presented in Figure[T2]to map specific exam- 




1 1.5 
Redshift 

Figure 12. The redshift sensitivity of the modified gravity parameter ■ 
The weight function (/"(z) is defined in equation j22| and evaluated with 
equation |23l. 



pies of theoretical models onto our parameter space, by evaluating 
equation {12\ . However as stressed earlier, we do not aim to place 
rigorous parameter constraints on any particular family of models. 



8.1 f{R) 

A more general form of the Einstein-Hilbert action replaces the 
Ricci scalar R with an arbitrary function f{R) such that 



S - 



//(«), 



-gd X, 



(24) 



where g is the determinant of the metric tensor. This defines the 
broad class of f{R) models. One of the most difficult tasks for 
any modified gravity model attempting to replace dark energy is 
to satisfy the stringent Solar System constraints, and most natu- 
ral choices of the function f{R) fail to do so. The subset of f{R) 
models which have attracted interest are those which employ the 
so-called chameleon mechanism, where departures from GR are 
strongly suppressed in regions where R is large, only emerging 
when R is sufficiently small. Our location within the potential well 
of the Sun and the Milky Way halo may be sufficient to shield us 
from this unusual gravitational behaviour. 

For a particular subset of f{R) models which are capable of 
satisfying Solar System tests, the departure from GR may be char- 
acterised as ( |Zhao et al.|2012[ > 



fi{k,a) 



3 + ?.{aM/kY 



(25) 



where the scalaron mass M(a) = 1/^3 d2//di?2. For any given 
redshift and wavenumber, the value of ^ lies in the range 
^ < |. This generically enhances growth, so we expect this 
family of models to lie vertically above the point (0, 0) in Fig- 
ure pi We parameterise M = MQa~'' and take as an example 
Afo ~ 0.02/iMpc~^ and cr = 3, corresponding to the type of 

(12012'). In /(i?) models the lensing 



Zhao et al. 



model explored in 
potential for a given mass distribution is unchanged from the case 
of GR, and SO ^ — 0. Our measure of is dominated by the 
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RSD data and by evaluating equation \22\ this model ought to map 
onto our (/x, E) plane at a value of /ig — 0.81. This is despite 
^(a) never rising above 1/3, but the stronger deviation from GR 
at high redshift compared to our parameterisation allows a higher 
effective value of /lo to be generated. Our data therefore disfavour 
the more aggressive f{R) models, in agreement with the findings 
of |Zhao et allpOlIl i. 

8.2 DGP 



The higher-dimensional model proposed by |Dvali, Gabadadze, &| 
|Porrati| ( |2000^ has proved one of the most popular gravitational 
models for theoretical and experimental exploration. It is of particu- 
lar interest due to its physical motivation and specific observational 
prediction, namely a value of the growth index 7 = if > or equiv- 
alently /xq ~ "H^a; Sq = 0. However in its most natural form 
there appears significant disparity with observational data (see for 
example |Lombriser et al.|2009 ). This can be bypassed by introduc- 
ing a cosmological constant, although this somewhat compromises 
the motivation of the DGP model. 

In general the behaviour of DGP is a scale-independent sup- 
pression of the growth of structure. Like f{R), the lensing potential 
is unaltered so Eq"^^ — 0. Adopting a fiducial value SIa = 0.27 
leaves us with /Xq"'^ — —0.38. This negative value of fi remains 
consistent with our data. 



8.3 Newtonian Limit 

The point at fio = 0; Eo = —0.5 corresponds to a scenario where 
the present-day cosmological gravity resembles Newtonian gravity. 
That is, the angular deflection of light is half the value predicted by 
Einstein at 2: = 0, while the forces experienced by galaxies remain 
unchanged. This model was famously rejected in favour of General 
Relativity by |Dyson et al.[ ( |192 0l. The difference here is that we are 
performing this test on cosmological scales. Whether we combine 
our gravitational data with either BAO or the Ho prior, the results 
are the same, a value of Eq — —0.5 is strongly disfavoured by the 
data. 

Another regime of interest is that of fio < — 1, for this corre- 
sponds to the anti-gravitational behaviour of galaxies, in that they 
are repulsed from clusters and attracted towards voids. One might 
naively expect that if dark energy is gravitational in origin, then the 
repulsive nature of the expansion may also be seen within the per- 
turbations. This region is disfavoured by the data when assuming a 
flat ACDM background, but this is no longer the case for more gen- 
eral backgrounds since the regime w < —1 permits more negative 
values of fio, as seen in the central panel of Figure|6] 



8.4 Interacting Dark Energy 

When conducting cosmological tests of gravity, in this work and in 



other s ([Daniel et al.|20l"0{|Bean & Tangmatitham|2010||Zhao et al 



20T0l|Song et al.|2011[|Reyes et al.|2010[|Tereno et al.|2011[|Zuntz 



et al.||20f2| " Zhao et al.||2012[ Rapetti et al.||2012| , we necessar- 
ily make certain assumptions such as the lack of anisotropic stress. 
One assumption not often made explicit is that dark energy exhibits 
negligible interactions with dark matter, aside from their mutual 
gravitational repulsion. If the two fluids were to exchange energy 
and momentum, this readily alters the growth of structure ( |Amen-| 
|dola,2Q00, Simpson et al.|2011[|Simpson|2010^ , and so could gen- 
erate a false signature of modified gravity. 



In the simplest interacting models, we would expect to see an 
alteration in the growth of large scale structure, but no direct change 
to the lensing potential. Thus detecting a non-zero value of Eq is a 
prerequisite for confirming a truly modified gravity scenario. That 
is to say, we cannot perform a cosmological test of gravity without 
gravitational lensing. 

As an illustrative example, we consider the scenario of dark 
energy decaying into a form of dark matter. If this results in the dark 



matter density following a power law p — p^a 



-3 + e 



the growth rate 



of cosmic structure is subject to a constant decrement ( [Simpson] 
|etal.|2011| l 



67 
55' 



(26) 



Large values of e lead to a significant quantity of dark energy in 
the early Universe, and as such are disfavoured by observations of 
the CMB. Taking e = 0.01, which is small enough not to perturb 
the primary anisotropics of the CMB, leads to an effective value 
^0 = —0.99 and Eg = 0, a region which is disfavored by our 
results. 



9 CONCLUSIONS 

In this work we find no indication of a departure from General Rel- 
ativity on cosmological scales, having explored the behaviour of 
non-relativistic and relativistic particles in the WiggleZ, 6dFGS, 
and CFHTLenS data. Both the motions of galaxies and the deflec- 
tion of their emitted light conform to the predictions of Einstein's 
theory, at the level of precision permitted by current data. We place 
limits on our gravitational parameters /io = 0.05 ± 0.25 and Eo = 
0.00 ± 0.14 for a flat ACDM expansion history. This corresponds 
to deviations in the present-day time dilation and spatial curvature 
of A*/*GR = 0.05±0.25and A<I>/$GR = -0.05 ±0.3 respec- 
tively. However these errors are significantly correlated, as shown 
in Figure [TT] 

For ground-based and solar system tests, where light-travel 
time is short, it is relatively straightforward to conduct a geometric 
measurement. Cosmological experiments do not share this luxury, 
and considerable effort is required in order to determine the spa- 
tial positioning of the test subjects. The precision at which we map 
the geometry of the Universe remains one of the limiting factors in 
determining its gravitational potentials. This is particularly appar- 
ent when increased freedom in the cosmic expansion history is al- 
lowed, arising from a variable dark energy equation of state or non- 
zero global curvature. Further geometric data could be included in 
the analysis, such as those derived from the light curves of Type la 
supernovae. However combining a large number of datasets leaves 
us increasingly exposed to the effects of systematic errors. For the 
purposes of gravitational lensing, another source of uncertainty is 
the form of the matter power spectrum on nonlinear scales, and this 
is an issue which will likely become more acute with the advent 
of surveys spanning a large fraction of the sky. Higher resolution 
N-body simulations can help to a certain extent, but modelling the 
influence of baryons may prove particularly problematic. 

The absence of a compelling theoretical rival to GR makes a 
suitable choice of gravitational parameters rather uncertain. In or- 
der to minimise the number of degrees of freedom, and to main- 
tain sensitivity to a broad range of models, we adopted a scale- 
independent modification to the Newtonian and curvature poten- 
tials. The scale-dependence of our results is implicit in the wide 
range of sensitivities each probe possesses. The cosmic shear signal 



© 0000 RAS, MNRAS 000, 000-000 



14 F. Simpson et al. 



from CFHTLenS is sourced by density perturbations of the order 
~5 /i~^Mpc, while the redshift space distortions from WiggleZ in- 
volve Fourier modes with wavelengths over 30 /i^^Mpc. Further- 
more the ISW effect causes the largest angular scales of the CMB to 
receive significant contributions from gravitational perturbations in 
excess of ~500 /i~^Mpc. If any inconsistency were seen between 
these results this could have been an indication of scale-dependent 
behaviour. Our reasons for focussing on a time dependent varia- 
tion that mimics the dominance of dark energy rather than using 
constant ^ and E values are twofold. First our motivation for seek- 
ing modification to General Relativity stems from dark energy, a 
phenomenon which only appears at late times. Furthermore, if we 
were to investigate constant fj, and E parameters, our constraints 
would be dominated by the primary anisotropics of the cosmic mi- 
crowave background, which as shown in previous studies ( |Bean| 
|& Tangmatitham|2010[ |Zuntz et al.|2012^ is consistent with a GR 
framework. 

Solar system experiments have already eliminated many po- 
tential modifications to GR which could have provided a natural 
explanation for the accelerating Universe. Here we have begun to 
disfavour some of those which remain. While the volume of param- 
eter space available to modified gravity models continues to shrink, 
the cosmological constant retains its position as the most viable so- 
lution to the dark energy problem. 
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APPENDIX A: TREATMENT OF THE CMB 

In our analysis we use the I ^ 100 WMAP7 data from the Cf ^ 
and Cj^ power spectra to constrain geometric cosmological pa- 
rameters, as well as a benchmark for the amplitude of matter per- 
turbations. In contrast, the I < 100 multipoles provide additional 
information on the gravitational potentials through the ISW. Dis- 
carding the low multipoles would severely compromise our ability 
to constrain the primordial amplitude of density perturbations. This 
is because on small angular scales there is a strong degeneracy be- 
tween the optical depth r and the amplitude Ag, since both param- 
eters induce a simple scaling of the amplitude. To circumvent this 
problem we first run a likelihood analysis on the low multipoles 
2 ^ ^ < 100, marginalising over our full parameter set includ- 
ing /lo and Eq. This establishes a prior on r, which we then use 
in our analysis of the data, except in those cases where the low 
multipoles are included. Despite the addition of modified gravity 
parameters, we find little degradation in the constraint on r com- 
pared to the standard WMAP7 result which assumes GR. This is 
due to the unique signature left by r at the large scales of the Cf^ 
spectrum, which cannot be replicated by the other parameters. 

To illustrate the impact of removing the WMAP7 data, we 
present in Figure [AT] the results for CFHTLenS +H0 in red, and 
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Figure Al. Constraints on our modified gravity parameters ^nd Eq from 
redshift space distortions (green), weak lensing (red), and combined (blue), 
for a flat ACDM background. A prior on Hq is included, but there is no 
additional information from the CMB. The dashed contours represent the 68 
and 95 per cent confidence intervals, while the cross denotes the prediction 
ofGR. 



WiggleZ +-ffo in green. The substantial degradation in these con- 
tours highhghts the importance of including the high— ^ data. 
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